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the shell. This program was used to determine the 
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nondimensional inplane stiffness 

the modulus of elasticity of the shell 
the rise of the spherical cap at the pole 
the thickness of the shell 

the mass density of the shell 


the meridional bending moment per unit length 


the Fourier index 
a nondimensional applied load 


the nondimensional critical pressure 

the classical buckling pressure of a complete 
sphere 

a column matrix containing the coefficients 


of the = term in the series expansion of the 
applied load 

the normal distance from the axis of revolution 
to the surface of the cap 


the normal distance from the axis to the cap 


in the base plane; the maximum value of [r 


the radii of curvature in the s and@ 


directions, respectively 

the meridional distance along the surface 
of the shell 

the nondimensional time 

the time 


a reference time 





<i 


MAX 


(n) 
W 


§t 


(0) 


the displacements in the s, © and J 
directions, respectively 

nondimensional series coefficients of U, V, W 
a nondinensSsional mneasure of the volume of 

the shell deformation 

the peak in the time history of the 


parameter V 
the displacement in the ¥F direction 


th 
in the nn harmonic 
the nondimensional time increment 


distance between stations 


the nondimensional parameter governing the 
Magnitude of the load applied in the asymmetric 
harmonics 

the coordinate normal to the surface of the shell 
the circumferential angle measured about 

the axis of revolution 

a nondimensional geometric parameter used 

to describe the spherical cap 

Pcisson's ratio 

the normal distance from the base plane to 


the middle surface of the undeformed cap 
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I. INTRODUCTION 


In 1973 a digital computer study was presented by Ball 
and Burt {1} for the dynamic buckling load of clamped 
Shallow spherical shells subjected to axisymmetric and 
nearly axisymmetric step-pressure loads. A static buckling 
analysis of the same spherical shells had been. carried out 
in 1970 by Stilwell and Ball [2]. In these two studies the 
digital computer program SATANS-I [3] was used to calculate 
the critical buckling pressures for a large range of shell 
sizes. Other studies of the puckling of shallow shells have 
been conducted by Huang {4,5], by Stephens and Fulton [6], 
by Lock et al. [7], by Stricklin [8], and most recently by 
Akkas (9). In Reference 1 the results from these other 
studies, except for those by Akkas, are compared with the 
results from SATANS-I for both static and dynamic buckling. 
In the axisymmetric static analysis tae comparison with the 
results obtained by Huang [4] revealed that the SATANS-I 
results were higher than Huang's results for several shell 
sizes. In the dynamic, axisymmetric buckling analysis the 
SATANS-I results again either agreed closely with, or were 
somewhat higher than , the results by Huang [5], Stephens 
and Fulton (6), and Stricklin [8]. However, it was noted 
then that there was a general lack of consistent agreement 
among any of the sets of results. As a consequence, it 
appeared at that time that the axisymmetric buckling problem 
had not yet been totally resolved and that additional 


studies would be appropriate. 


In the asymmetric dynamic buckling analysis of Reference 
1 the few comparisons that could be made for the critical 
load also indicated that the SATANS-I results may be too 





high. A ccmparison of the recent estimates for the 
asymmetric dynamic buckling load obtained by Akkas [9] with 
the SATANS-I results also reveals the SATANS-I results to be 
well above those of Akkas [9]. However, it should be noted 
that the results obtained by Akkas were from his attempt to 
obtain a lower bound on the critical asymmetric load. This 
bound on the buckling load 1s obtained without the execution 
of a complete transient response analysis on the asymmetric 
part of the response of the shell, as is done in SATANS-I. 
In Akkas* analysis (Problem 1) the transient nonlinear 
axisymmetric response is computed, and a determinant is 
examined for possible bifurcation into asymmetric moticn at 
each time step. The minimum load at which the determinant 
becomes zero is defined as the lower bound of the critical 
load. 


As a consequence of the generally high buckling loads 
predicted by SATANS-I, a re-examination of the static and 
dynamic buckling of the shallow spherical shell was made in 
an attempt to determine the possible cause, or causes, of 
the high buckling loads. In our search we discovered that a 
modification of the manner in which the pole conditions are 
numerically approximated significantly lowered the buckling 
loads to values that are now in good agreement with the 
other results. The new procedure for handiing the pole 
condition is given in section III of this thesis. The new 


buckling results are given in section V. 


In addition to the pole condition modifications and the 
new buckling results the author has also made another 
Significant change to the SATANS family of codes. [In 
particular, the SATANS-II program for the geometrically 
nonlinear analysis of totally arbitrarily loaded shells of 
revolution, developed by Ryan [10] in 1972 to handle more 
complex and larger problems, was modified to make the 


computer memory reguirement a variable quantity. This 
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quantity is specified by the user to fit the particular 
problem being run. It eliminates the large core requirement 
of SATANS-II for small problems and allows for much larger 
problems to be solved than could be solved by SATANS-II. 
The new program with the pole condition and memory 
modifications will hereafter be called SATANS-IIA. ines 


described in section II. 
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SATANS-II was developed by Ryan [10] from SATANS-I and 
incorporated the full trigonometric expansion of the applied 
load and solution vector, and introduced the handling of 
imperfections into the code. These modifications allow the 
analysis of shells under totally arbitrary loads, as well as 
imperfection studies on actual shells with measured 
imperfections {11]. Unfortunately, the original deck of 
cards for SATANS-II was destroyed. Professor Johann Arbocz 
of CALTECH had a listing of SATANS-II and punched a deck of 
cards with the changes to SATANS-I given in that listing. A 
copy of this deck was sent to Professor Ball. These cards 
have been added by the author to the original SATANS-I 
described by Ryan {10} and a complete version of SATANS-II 
has been reconstructed. SATANS-IIA is a modification by the 
author of the reconstructed SATANS-II program. A listing of 
SATANS-IIA can be found in Appendix A. The listing contains 
an example problem for the dynamic analysis of a clamped, 
truncated cone subjected to an impulsive loading which is 
uniform along the meridian and varies en a cosine 
distribution over one-half orf the circumference. This 
problem is a sample problem suggested by the Lockeed 
Missiles and Space Corp. [12]. A condensed version of the 
output from the example problem iS given in Appendix B. 
Input data preparation for SATANS-IIA can be found in 
Appendix C. The basic uSers manual, which includes 
preparation of input subroutines and the theory of the 
program, is contained in Reference 3, which can be obtained 
through COSMIC (M70-10098, LAR-10736), or ASIAC [13]. A 
users manual which includes preparation and handling of 


imperfection data within the SATANS programs can be found in 
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Ref. {10}. The above information, along with the following 
discussion, will inform the user on the capabilities and 


proper use of SATANS-IIA. 


The modification of the SATANS-II program to make its 
core requirement variable was accomplished by putting in a 
Single dimension statement at the beginning of the progran, 
with subsequent dimensioning within the subroutines to only 
the first element of the vector or matrix. This is a 
convenient feature of the FORTRAN-IV language in which the 
program is written. The actual vector and matrix sizes are 
transmitted to the subroutines by. an individual parameter 
list. Construction of the initial dimension statement and 


core request size is as follows: 


The basic size of the program on the IBM-360/67 Digital 
Computer, without the initial dimension statement, 1s 
272,000 bytes. This figure includes approximately 19,000 
bytes of buffer space required for execution. Within the 
Main dimension statement are fifteen variables. However, 
only three parameters are needed to specify the sizes of 


these fifteen variables. 


Let a= The number of stations along the meridian of the 
shell times the number of harmonics considered. 


Let b= a, plus two fictitious stations times the number 


of harmonics considered. 
Let c 


The number of harmonics considered. 


The main dimension statement would then be constructed 
as, 
DIMENSION P(4,4,a), DEE(4,4,a), DST(4,4,a), X(4,a), 
PHIXB(a), PHITB(a), Z(4,b), ZO(4,b), 
[eo Lo ae), ZDOT (4,b), 1S(99,c), 
Someone, koe s,c) , JD(99,c) 
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The 99's above limit the user to 99 harmonics in any one 
run and an unlimited number of meridional stations. The 


core requirement for the general case would be, 


272,000 + 216a + 80b + 1584c = bytes of core reguired. 


For a sample calculation of the core reguirements 
consider the example of a spherical cap with 40 stations 
along the meridian, and an asymmetric analysis with two 
harmonics. MTherefore, 


a= 4O(stations) x 2(harmonics) = 80 
b= 80 + 2 x 2(harmonics) = 84 


c= 2 (harmonics) 


Thus, for the variables P, DEE, DST, 
3 x (4x4x80) = 3840 (words) x 4 = 15,360 bytes 


for the variable X, 
4 x (80) = 320 (words) x 4 


1280 bytes 


for the varibles PHIXB, PHITB, 
2 x (80) = 160 (words) x 4 = 640 bytes 


for the varibles Z, ZO, Z2, 23, ZDOT, 
5 x (4x84) = 1680 (words) x 4 = 6720 bytes 


mastiy, for the varibles ID, JD, IS, JS, 
G4 x (99x2) = 792 (words) x 4 = 3168 Eytes 


Therefore, the total size of the main dimension 
statement would be 27,168 bytes. This figure would be 
rounded up to the nearest even thousand bytes, ie. 28,000 
bytes. Finally, the core reyguirement for this example 
problem would be 


272,000 + 28,000 = 300,000 bytes. 
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Tift. IMPROVED POLE ROUTINE 


= ee ee a = aes £= ae 2 EE SS 


The SATANS code is based upon Sander's geometrically 
nonlinear equations under the conditions of small strains 
and moderately small rotations. The formulation is in four 
second order nonlinear partial differential equations in 


terms of U, V, W, and M, where U, V, and W are the 
S 


meridional, circumferential and normal displacements 

respectively, and M is the meridional bending moment. The 
S 

nonlinear partial differential equations in the coordinates 


S, 9 , and t are reduced to uncoupled sets of linear 
differential equations in s and t by expanding the variables 
in trigonometric series in the circumferential coordinate @, 
and treating the nonlinear terms as pseudo loads. The first 
and second derivatives in the meridional coordinate s are 
replaced by the conventional central finite difference 


approximations, ie. 


ye 1728 (fz) = {Zz}. .) (1) 
1 1+1 1-1 
and 
. 2 
{2}. = 1/0 ({Z}. eee je {Z}_ 2) (2) 
1 1+1 1 1-1, 
th 
where {z}. is the vector of JU, V, W, and M at the i 
1 Ss 


station, 4 is the uniform dimension between stations, and 


primes denote partial derivatives with respect to s. 


Applying these approximations to the governing set of domain 
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eguations leads to 


(Cj. {Zz}. Capo ietzie et fA} {2} . Seo ( 3) 
ae 1-1 1 1 a 1+1 al 


When the shell does not have a pole, fictitious stations 
one increment off of the shell are introduced at each end. 
Both the governing domain equations and the boundary 
conditions are applied at the two boundary points. Thus, 
all finite difference approximations to the derivatives, 


including those of the boundary conditions, are of order 
2 
a . However, prior to the development of SATANS-iita, the 


treatment of the conditions to be applied at a poie at 
either end of a shell was handled by a simple Euierc rorward 
or backward difference approximation to the PAcse 
derivative, with truncation error of order 4. For example, 
for a pole at s= 0, where i= 1, the first derivative at the 


pole waS approximated with 


=e ~ ; (4 
(z¥ /O& ({z}, - (23) ) 


At the time this procedure for handling the pole 
conditions was developed (1967) it was thought that this 
would not significantly alter the solution. However, it has 


since been discovered that such 1s not the case. 


For the new pole routine, an expanded forward difference 


2 
approximation cf order A is used at s= 0 which takes into 
account the two stations after the pole, instead of just one 
Station after the pole as in the Euler schene. This 


approximation is 


{z} = 1/24 (-3{z} + 4{z}_ - {z}). (5) 
1 1 2 3 
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The 


variables at a pole are derived in Reference 14. 


For N= QO, 


Applying equation (5), these conditions can be put 


matrix form 


100 0 u 
010 Q Vv + 
9 0-3 0 Ww 
mo 0-3 m 

s’ 1 


where the above 3 matrices 


SATANS programs. 


For N= 1, 


where the plus 


minus sign at a final pole. 


Songditions is 


-3 00 0 
141 0 0 
00 10 
000 1 


For N=2, 


the matrix form is 


0 0 


lo Oo OO 
B « 4 £& 


FOr 


conditions 


u 


=V 
1 


to 


1 


pe 


é é 


=w° =m 
1 


s1 


000 0 u 
00 0 Q V 
004 0 wW 
000 4 m 
are DL 
u+#v=2s=uwew 
1 1 
sign applies 
Th 
uQdQd dQ 
+ 000 0 
000 0 
000 0 
1 
u= v= w= m = 0 
s 
000 0 u 
000 0 V 
000 0 W 
000 4 m 
u= V= w= mn = 
s 


7 


imposed 


QO. 


S 


a 


at 


e 


B x << ¢ 


upon the dependent 


They are 3; 


into the 
000 0 u 
+ 00 0 0 Vv =i" 
0 0-1 0 W 
00 0-1 m 
zZ S3#+3 
DG, and DF within the 
n = QO, 
s 


an initial pole, and the 


Matrix fcrm for these 
a 0 1] 

+ 0 0 Vv = 0. 
0 0 Ww 

2 ae I a 3 
00 0 0 u 

+ 00 0 0 Vv = Q. 
00 0 0 wy 
00 0-1 m 

Z s° 3 





and DL= identity matrix, DG= DF= null matrices. 


The sclution procedure in SATANS is an eiimination 


scheme and starts with 


es = Pd, (23, . ie ’ (6) 


where the values in cP], based upon the Euler approximation 
are defined in Reference 14. The higher order approximation 


defines a new ao This new ee 1s obtained by 


Simultaneously solving the pole conditions 
[DL] (23. + [DG] at + [DF] Soe = {0}, (7) 


and the domain eguation at station 2 next to the pole 


C B A = ; 8 
[ se va tz + [ 1), {3}, (8) 
to eliminate Thus, 
ane = EG - {B 2 9 
a J, ({93 , [ J, (2} [ J, {2} )) (9) 


Substituting equation (9) into equation (7) gives 
| 
DL DG DF A =e = 
f be oot zit ee dee eee 


B = 0. 1 
[ a (23) (10) 


Combining like coefficients of the {z} vector leads to 
Sel eeonrA Le) ) 2) + ( [Dc] - (DF) fay fa) ) 
2 2 1 2 2 


eee A). : 11) 
(23, euler) [ J, {9}, ( 


Finally, solving for (2) yields 
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—_ 


Zz fDi =hre rah x °C DG - DF x A! x B Zz 
{ i [ 5 =! [ ; ae { Ne 


L- DP x A’ xc =" Dex uA : 12 
+ { D eS Senos (12) 


=1 
Thus P =- (DL - DF x A“? x C DG - DF x A’ xB and 
PL J, [ 5 a [ 5 5! 


-1 
x = [DL ~- DF x A’ xc - DF x a’! . The new 
{x}, [ 5 ae [ 54 (g}, 


{ Et, matrix has been placed into the "PMATRX" subroutine of 
SATANS-IIA and the new at vector has been placed in the 


"FORCE" subroutine. 


A listing of the pole routine may be found in Appendix 
D. To inccrporate this new routine into a SATANS-I or-II 
program, first proceed to the "PMATRX" subroutine and remove 
the fifteen cards that are between, but not including, 
“TF(NN.GT.2) GO TO 90" and "11 CONTINUE". These cardS are 
located after statement number "14" and just before 
statement number "11". Replace the cards removed by the 
ones listed in Appendix D which read from" C IN PMATRX " to 
“ 90 M3=MN". Then proceed to the "FORCE" subroutine and 
remove statement number "10". Replace statement number "10" 
with the nine cards listed in Appendix D which read from " C 
mee rORCE™’ to "DO 11 I= 1,4" . Also piace “ COMMON 
PEeho/EBCINL, IBCFNL " into the common area of the " FORCE " 


subroutine. 


This completes the implementation of the new pole 
routine into either SATANS-I or II. 
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IV. PROBLEM DESCRIPTION 


= ap SS Sw ae ae a ee 


The geometry of the shallow spherical shell used in this 
study is identical to that used in Reference 1. BEwerley; 
the shallow shell can be specified by the non-dimensional 
parameter X , where 

2 1/74 1f2 

A = 2 3 ( ve a) ( H/ hb) : (1) 
H is the rise of the shell, h is the thickness, and Y is 
Poisson's ratio. The mass density of the shell is ao. All 


Shells analyzed had the following dimensions; 


” 


Radii of Curvature R= Ro = 250 inches 
s 

Thickness h = 0.25 inches 

Modulus of Elasticity E = 30,000,000 psi 

Poisson's Ratio Y= 0.3 


All buckling pressures obtained will be listed as a 
percent of the classical buckling pressure of a complete 


sphere, ae where 


2 2 i772 
qa RS ae nce Pes t= Vv ) 7 (2) 
Forty stations were used over the meridian. The 
nondimensional time increment § t, where 
‘ 1/2 
Sey ae / &) ’ (3) 
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was taken as 0.05 for 3000 time steps, which 1s a total 
nondimensional time of 150. In addition, the axisymmetric 
analysis was repeated with a larger time step of 5 t= 0.2 
for a total time of 600. In this study m was selected such 
that t is equal to fT. The necessity for the long response 


time 1s explained in Reference 6. 


In the axisymmetric analysis only the N= 0 harmonic is 
considered. However, in the asymmetric analysis a_ second 
harmonic is excited by applying an incremental load in that 
harmonic. In addition, analyses of the shells X = 6, 7.5, 
and 11 were made using five harmonics. The step pressure 
load for the axisymmetric harmonic is 

(9) 

cies 3} = 2 a {1}, (4) 
and the step pressure load for the asymmetric second 
harmonic is 

(n) (n) 


ide j} = P Ty E {1}, (3) 


where n> QO, and ee is taken as 0.0001. The value taken 
for the second harmonic in the asymmetric analysis was the 
same as the critical naarmonic for the static buckling 
analysis presented by Stilwell and Ball [2]. When there was 
an uncertainty as to which was the critical static harmonic 
the two harmonics in question were both tested. Run times 
using SATANS-IIA with a two-harmonic analysis for 3000 time 
steps and 40 stations on the meridian took an average of 28 
Minutes on the IBM 360/67. 


The parameter used to determine the minimum load at 


which dynamic buckling occurs is the peak value of V, called 
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V , where V is defined as 
MAX 


"5 (0) fh 
ve rs ary | r¥ ar (6) 
0 9 


ris the normal distance from the axis to the shell, - 1s 
(9), . 
the maximum value of r, W is the normal displacement of 


the axisymmetric response and § is the vertical distance 


from the tase plane to the undeformed shell. The Visa 
measure or the volume or the shell deformation. The Fortran 


statements computing V and ee are given in Appendix E. 


When working a problem that requires these calculations the 


Mineteen cards are inserted directly into the "DYNAMIC* 

subroutine right after the "IF" statement that calis tke 
“OUTPUT” subroutine. 

For convenience, the response in each asymmetric 

0) 


harmonic is also measured uSing equation (6), with W 


(n) = . 
replaced with W - The parameter V for the asymmetric 


harmonics does not represent a volume of deformation as it 


does for the axisymmetric harmonic. It can, however, be 
used to indicate the relative excitation of the asymmetric 


harmonics. 


The buckling criterion for both the axisymmetric and the 
asymmetric dynamic buckling analysis defines the critical 


load as that load P where a very small increase in P causes 


‘a very large increase in V - This is the same criterion 
MAX 


ae, 





as that used in Ref. 


{ 1]. 
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Ve. RESULTS AND DISCUSSION 


= eet a Se = SZ 2 ea ee eee 


A. STATIC AXISYMMETRIC BUCKLING ANALYSIS 


Table I presents the new results from the static 
axisymmetric buckling analyses for A = 4 through 13 using 
the new pole routine. The two upper curves in Figure 1 
present a comparison of the new results obtained by 
SATANS-IIA with those obtained by Stilwell and Ball [2] 
uSing the SATANS-I program. As can be seen in this figure, 
fairly significant changes in the buckling load occurred in 
the neighborhood of X» = 4,5, and 9; and somewhat smaller 
differences cccurred in the region A = 10 through 13. The 
upper data points in Figure 2 present the comparison of the 
hew results from SATANS-IIA with those obtained by Huang 
[4 }j. This comparison shows a very good agreement between 
the two sets of results, except for the largest values of dX. 
The new results have eliminated the differences that existed 


between the SATANS-I results and Huang's results. 


B. DYNAMIC AXISYMMETRIC BUCKLING ANALYSIS 


Figure 3 presents the new results for the peak value of 


V versus P for the various values of A tested. Table 
Mea X 


If presents all of the new results for the dynamic 


axisymmetric buckling load. These loads are selected from 


figures constructed just like Figure 3. In every case, 
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except for »* = 4, a value of P slightly above the Pan : 
ir 


value caused a V me indicative of buckling, as well as a 
M 


nonconvergence of the iterative solution procedure. 


The lower two curves of Figure 1 present a comparison 
versus » of the new axisymmetric dynamic buckling results 
with the previous buckling results obtained by Ball and Burt 
{1}. In every case the new critical pressure is lower than 
the critical pressure obtained using the Euler approximation 
at the pole. 


The lower data points of Figure 2 present a comparison 
of the new results with those obtained by Huang [5], by 
Stephens and Fulton [6], and by Stricklin [8]. Just as in 
the case of the static axisymmetric buckling analysis, the 
new results compare much more favorably with the other 
results than did the results of Reference 1. Its 
interesting to note that the new results now tend to be 
Slightly lower than the other results, whereas the resuits 


of Reference 1 were higher for almost all values of A . 
Ca DYNAMIC ASYMMETRIC BUCKLING ANALYSIS 


Table III presents the new results for the critical 
pressures obtained from the dynamic asymmetric analysis. 
The second harmonics, or critical static harmonics, used in 
the analyses are also presented in Table III. A comparison 
of the critical pressures from the asymmetric analyses, 
Table III, with tne critical pressures from the axisymmetric 
analyses, Table II, reveals that only the shell A =6 
buckled at a load below the axisymmetric buckling load. For 
the shell A = 7 the critical buckling load was slightiy 
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larger when asymmetric motion was considered. In all other 
cases the buckling was not influenced by the presence of the 
second harmonic. These new buckling results and those by 
Ball and Burt [2] are plotted in Figure 4&. The new resuits 
can be seen to be significantiy different from the SATANS-I 
results, where the asymmetric buckling loads were lower than 
the axisymmetric loads for five out of the ten values of 
tested. 


Except for » = 6 and 7, the relationship between ae 


and P for the N= 0 harmonic, in the two-harmonic analyses, 


was found to be essentially identical to the relationship 
found in the axisymmetric buckling analysis shown in Figure 


ee Table IV A presents the V versus P data for both the 
MAX 


N= 0 harmonic and the second harmonic, for all values of A 


tested, except for A» = 6. Note that, except for X#=7 and 


a, Vv for the asymmetric harmonic is generally very 
M 


small, even when the V¥ for the N= 0 harmonic indicates 
MAX 
that the shell has buckled. Thus, except for the shells A= 


6 and 7, the presence of the asymmetric motion does not 
influence the axisymmetric motion, and except for the shells 
» = 6, 7 and 11 the asymmetric motion is very small prior 


to buckling in the axisymmetric harmonic. 


A more detailed analysis of the shell X = 6 has been 
conducted since it was the only shell that revealed any 
Significant axisymmetric sensitivity to asymmetric motion. 
This shell was studied using two two-harmonic analyses (N= 
O, 1 and N= 0, 2) anda five-harmonic analysis (N= 0, 1, 2, 
3, and 4). Figure 5 and Tables IV B and IV C contain values 


of V versus P for both of the asymmetric harmonics, N= 1 
MAX 
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and N= 2, in the two two-harmonic analyses, as well as the 

values of ns for the axisymmetric harmonic, N= 0. Figure 

6 and Table IV D present the values of ve versus P for the 
iA 


N= 0,1,2,3, and 4 harmonics from the five-harmonic study. A 
comparison of the critical buckling load predicted from the 
results of the two two-harmonic analyses in Figure 5 with 
the critical load from the five-harmonic analysis obtained 
from Figure 6 shows that the presence of the additional 
harmonics results in the shell buckling at a slightly lowec 
load (0.50), with significant motion in the N= 1 harmonic 
instead of the N= 2 harmonic (see the norconverged solution 
at P= 0.51), which is the critical harmonic for static 
asymmetric buckling. Studies using five harmonics have also 
been conducted for A = 7.5 and A= 11. As can be seen in 
Table IV D the critical harmonic for A= 7.5 remained N= 3; 
however, significant motion occurred in that harmonic at P= 
wit and .44. In the case -of A = 11, relatively large 
asymmetric motion occurred in the asymmetric mode of N= 5 


vice 6 at a value of P= .46. 


The comparison of the new results for the critical 
pressure for dynamic asymmetric buckling with those obtained 
analytically by Stexek lan eer, by Akkas [9], and 
experimentally by Lock et al [7] is iliustrated in Figure 7. 
The comparison reveals an agreement with Stricklin in every 


case, in general a higher value of Beng than those obtained 
T 


by Akkas, and most importantly a very good agreement with 


Lock's experimental results. 


When making the comparison between the new results and 
those obtained by Akkas, it is necessary to look at _ the 
ditferences in the problem solution parameters used in the 
two studies. For example, buckling results obtained fron 


SATANS-IIA using the same time increment as used by Akkas, 
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§ t= .2 for 3000 time steps, were significantly higher 
than those using the time step of St= .05 for many values 
of A. Furthermore, the new results had , in some cases, 
instances of buckling occurring as far out in time as 130. 
Akkas, to shorten computer run times, observed the cap only 
for a time of less than 5. Purthermore, only the harmonics 
femi or 2 or 3 were studied by Akkas for shells A= § 
through 12. If the critical harmonic is not studied, the 
predicted load will be too high. Thus, it appears that 


Akkas' lower bound loads may not be true lower bounds. 


Two additional features of the shell response should be 
noted. First, shells A = 6, 7.5, and 11 exhibited a 
non~buckled response in the axisymmetric harmonic to a_ load 
larger than the defined critical buckling load. This can be 
seen in Tables IV A and IV C. Second, and most importantly, 
the buckling load proposed by Ball and Burt [1], and used 


here, defines buckling to occur when the V in the 
MAX 
axisymmetric harmonic undergoes a large change due toa 


small change in P. Another criterion for dynamic buckling 
in the asymmetric analysis discussed in Reference 1 is to 
define the buckling load as that threshold load that 
initiates significant growth in the asymmetric harmonic. 


Re-examination of the ¥V versus P data in Table IV A 
MAX 


through CD reveals that shells A= 6, 7, and 11 exhibited 


relatively large asymmetric motion at loads smaller that the 


defined buckling load when compared with other V values 
MAX 


for those shells, even though the numbers themselves were 
\ 


small when compared with the axisymmetric harmonic. Shellsa 
77.5 and 12 appear to be borderline cases. If the alternate 
criterion for buckling is used, the critical buckling loads 
for shells A= 6, 7, and 11 become 0.47, 0.45, and 0.45, 
respectively. The shells » = 7.5 and 12 could have buckling 
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loads as low as 0.40 and 0.44, respectively. These values 
are more conservative than the definition based upon 
axisymmetric response. These five shells are the same five 
shells that exhibited an asymmetric buckling load lower than 


the axisymmetric buckling load in Reference 1. 
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VI. SUMMARY AND CONCLUSIONS 


A digital computer program for the geometrically 
nonlinear. analysis of totally arbitrarily loaded shells of 
revolution (SATANS-II) was modified to more accurately 
account for the conditions at the pole of the shell. This 
program, called SATANS-IIA, was used to determine the 
buckling load of shallow spherical shells of various’ sizes 
when subjected to static axisymmetric, dynamic axisymmetric, 
and dynamic nearly axisymmetric step~pressure loads of 
infinite duration. The cap sizes ranged from >= 4 to 13 
including » = 7.5. A comparison was made between the new 
buckling results with the improved pole handling routine and 
the results that did not have the new pole routine. The 
comparison revealed a Significant change in buckling 
pressures, due solely to the change from an order A _ finite 


difference approximation of the first derivatives at the 


2 
pole to an approximation of order A . These new critical 


pressures are in very good agreement with the results fron 


other studies of the same spherical shells. This good 
agreement with other results, which came about as a result 
of the modification of the pole handling routine, is a 
strong aindicatzon that the manner ‘in which the pole 
condition is handled is vital to the accuracy of the 


solutions obtained. 


In the asymmetric analysis, two harmonics were included 
for most of the shells; the axisymmetric harmonic and one 
asymmetric harmonic. Five-harmonic analyses were conducted 


for three of the _  shelis. Two buckling criteria for the 
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asymmetric aralysis were considered. One defined buckling 


as that threshold load that caused a large increase ina 


deformation farameter, Vee in the axisymmetric harmonic. 


The other, more conservative than the first, defined 
buckling as that threshold load that caused a large increase 


in the ve value for the aSymmetric harmonic. Both values 
X 


have been presented. 


The new static axisymmetric, dynamic axisymmetric, and 
even the dynamic asymmetric critical buckling pressure loads 
appear to te fairly reliable results for perfect, shallow 
shells. The effect of realistic imperfections remains to be 


determined. 
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Figure 5 - PEAK DEFLECTION VERSUS P FOR THE ASYMMETRIC 
ANALYSES OF A= 6 (N=0,1 AND N=0,2) 
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A. TABLES 








3. FABLE fil Critical step-pressure ioads from the 
dynamic asymmetric analyses and critical asymmetric 
harmonics. 
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1. TABLE IV A. Two-harmonic analyses for all values of 





mee= 5 N= 0 and 2 N= 0 and 1 
a a 







N= 0 |.1659 | «1676 | .6606 |] w= 0] «1675 [.6606_ 
N= 2 |.0004787 | .0000566 | .0687 || n= 1|.0003145].0687_ 
Ree ee 


[oo] 7683 
[a= 1.001082 
= 7, N= 0 and 3 


psa ee 
n= of. 09452 [.09571 [09812 | .1005 | .1029| .1052 | .1099 
prof aaa [ae [3 


= N= 0 and 3 


a EC CY 
IN= 0|.0703 .07228| .07429 .2636| .07837 . 2076 
N= 3|.0001094|.001296|.0004304 | .002754 | .001188 | .000338 


- 46 
N= o| .200 




























4Q 





A = 8, N= 0 and 4 


; 


A = 9 N= 0 and 4 N= 0 and 5 


40 Je | 40 
.04738 |.04875 1576 | w= 0] .05012 
.00003597|.00004635|.00004497 || N= 5|.00008385 


10, N= O and 5 







= 
il 






















N= 0 
N= 4 


i = 
N= 0} .03239| .03347{ .1086 |.1217 .1288 | .1235 
N= 5]|.0000281 |.000 0472 |.00004125|.00002103|.0000449|.000114 


fee= 11, N= 0 and 6 


ps | ue ee 
N= o| .03910 | .04004 | .04099 | .09814 | .o4241 | .08824 
n= 6| .004595 | .01332| .02232| .02864 | .03955| .02813 


mee= 12, N= 0 and 7 





















N= 0/ .03236 | .03316 | .08633 
.00004214 | .0004561 | .00005158 


he = 13, N= 0 and 8 












I= 0} .02119 | .02185 . 06637 07381 
n= 8{.00001148].00001134 |.000006607].000008245|.000119 
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2. TABLE IV B. Two-harmonic analyses with N= 0 and 1, 
A = 6 cnly. 


[we | es] 0d —idzS Cd 
r= o| 1350 | .1305 | swat |. teeo | 1099 | .0us3_ 
T= 1],0002797].000195| .000205| .000926| .ouoet| «3068 


3. TAELE IV C. Two-harmonic analyses with N= 0 and 2, 
rN = 6 only. 


[swe .a5 | se] a7 | we | a9 
sof vanaf vaso fave [-isz0 [isso i305 | 79) 
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GY. TABLE IV D fFive-harmonic analyses for selected 


















snells. 
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p [7 [we | [so 
pa 01.133__[o1947 [1392 _| tuo _{-3797_{ tase) 2200 
= 1}.00021 |.02108 |.02215 | .003676|.3743_|.01276 |.0616_ 
Sibeenior—fasoss |e -asor _fersez [sate |. zn05 
| 00163 | 
os 
N= 0].3854 
N= 1|.3953 
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P75 °N= 0, 1,2,3,4 and 4 


Orgs S07 220 07429 2592 07837 2544 


00004855) .00004198 |.00006093 |}.0002737 }.0001167) .005952 







= 
“i 


2).0001164 00007456 |.0004184 |.0000982 }/.000788 |.0003188 
3 }.0001277 | .001187 0004597 |.0002853 {.00107 |.0003188 
41.0008224 | .0001898 ;.0002448 |.0000526 |.000280 1000134 


pee= 11 N= 0,4,5,6, and 7 


Pp [.s [us far [we [wa | 50 
w= 0] .03910 |.ouoo Lousa  Lowi95 | .ow291 | 1040 
i= 6|.002560 |.007020 | .02330 [02767 | .02602 | .o264u 
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APPENDIX A 


LISTING OF SATANS-IIA 
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